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We determine the sensitivity to a possible variation of the proton-to-electron mass ratio \x for 
torsion- wagging-rotation transitions in the ground state of methylamine (CH 3 NH 2 ). Our calcu- 
lation uses an effective Hamiltonian based on a high-barrier tunneling formalism combined with 
extended-group ideas. The /^-dependence of the molecular parameters that are used in this model 
are derived and the most important ones of these are validated using the spectroscopic data of 
different isotopologues of methylamine. We find a significant enhancement of the sensitivity co- 
efficients due to energy cancellations between internal rotational, overall rotational and inversion 
energy splittings. The sensitivity coefficients of the different transitions range from —19 to +24. The 
sensitivity coefficients of the 78.135, 79.008, and 89.956 GHz transitions that were recently observed 
in the disk of a z — 0.89 spiral galaxy located in front of the quasar PKS 1830-211 [S. Muller et al. 
Astron. Astrophys. 535, A103 (2011)] were calculated to be —0.87 for the first two and —1.4 for 
the third transition, respectively. From these transitions a preliminary upper limit for a variation 
of the proton to electron mass ratio of A/i/fi < 9 x 10~ 6 is deduced. 

PACS numbers: 06.20.Jr, 33.15.-e, 98.80.-k 



I. INTRODUCTION 

Recently, it was shown that transitions between acci- 
dently degenerate levels that correspond to different mo- 
tional states in polyatomic molecules are very sensitive 
to a possible variation of the proton-to-electron mass ra- 
tio, \i = m p /m e . Kozlov et al. [1 showed that transi- 
tions that convert rotational motion into inversion mo- 
tion, and vice versa, in the different isotopologues of hy- 
dronium (H 3 + ) have coefficients ranging from —219 
to +11 [2 . Similarly, Jansen et al. 013] and Levshakov et 
al. [5] showed that transitions that convert internal ro- 
tation into overall rotation in the different isotopologues 
of methanol have coefficients ranging from —88 to 
+330. Here, the sensitivy coefficient, K^, is defined by 
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For comparison, pure rotational transitions have = 
— 1, while pure vibrational transitions have = — ~ and 
pure electronic transitions have = 0. 

Accidental degeneracies between different motional 
states in polyatomic molecules are likely to occur if the 
energies associated with the different types of motions 
are similar. In this paper, we present a calculation of 
the sensitivity coefficients for microwave transitions in 
methylamine (CH 3 NH 2 ). Methylamine is an interesting 
molecule for several reasons: (i) it displays two large am- 



plitude motions; hindered internal rotation of the methyl 
(CH 3 ) group with respect to the amino group (NH 2 ), and 
tunneling associated with wagging of the amino group. 
The coupling between the internal rotation and over- 
all rotation in methylamine is rather strong resulting 
in a strong dependence of the torsional energies on the 
K quantum number, which is favorable for obtaining 
large enhancements of the if M -coefficients 0] . (ii) Methy- 
lamine is a relatively small and stable molecule that is 
abundantly present in our galaxy and easy to work with 
in the laboratory. Recently it was also detected in the 
disk of a high redshift (z = 0.89) spiral galaxy located in 
front of the quasar PKS 1830-2110. 

This paper is organized as follows. In Section[nJ we in- 
troduce the effective Hamiltonian used for calculating the 
level energies in the vibrational ground state of methy- 
lamine. In Section [TTTI we derive how the constants that 
appear in this Hamiltonian scale with /i. Finally, in Sec- 



tion IV we use the Hamiltonian and the scaling relations 
to determine the sensitivity coefficients of selected tran- 
sitions. 



II. HAMILTONIAN AND ENERGY LEVEL 
STRUCTURE 

Methylamine, schematically depicted in Fig. [I] is a 
representative of molecules exhibiting two coupled large- 
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Figure 1. Schematic representation of methylamine and vari- 
ation of the potential energy of methylamine as function of 
the relative rotation, 7, of the CH 3 -group with respect to the 
amine group about the CN bond and the angle, r, of the two 
hydrogen atoms of the NH 2 -group with respect to the CN 
bond. The two large amplitude motions, corresponding to 
inversion (h,2v) and hindered rotation (/i3„) are schematically 
indicated by the arrows. Note that inversion of the NH 2 - 
group is accompanied by a n/3 rotation about the CN bond 
of the CH 3 -group with respect to the amine group. 



amplitude motions, the torsional motion of a methyl 
group and the wagging (or inversion) motion of an amine 
group. A combination of intermediate heights of the po- 
tential barriers with a leading role of the light hydrogen 
atoms in the large-amplitude motions results in relatively 
large tunneling splittings even in the ground vibrational 
state. On the right-hand side of Fig. [T] a contour plot 
of the potential energy is shown with the relative angle 
between the methyl and the amino group, 7, on the hor- 
izontal axis and the angle between the NH 2 -plane and 
the CN-bond, r, on the vertical axis. The methyl tor- 
sion motion is indicated with the arrow labeled by h 3v 
whereas the amino wagging motion is indicated with the 
arrow labeled by h 2v . From the contour plot, it is seen 
that amino wagging motion of the NH 2 -group is accom- 
panied by a 7r/3 rotation of the CH 3 -group about the 
CN bond with respect to the NH 2 group. Consequently, 
the amino wagging motion is strongly coupled to the hin- 
dered methyl top internal rotation resulting in a rather 
complicated computational problem. 

In Fig. [2] the lowest rotational levels of the ground 
vibrational state of CH 3 NH 2 are shown. The level or- 
dering resembles that of a near-prolate asymmetric top 
molecule. In addition to the usual asymmetric splitting, 
every J, K level is split due to the different tunneling 
motions. The internal rotation tunneling splits each ro- 
tational level into one doubly degenerate and one non- 
degenerate sublevel. Each of these sublevels are further 
split into two due to the inversion motion. Together, 
this results in eight levels with overall symmetry Al, A2, 
Bl, B2, El + 1, E2 + 1, El - 1 and E2 - 1 for K > 
and four levels for K = 0. The +1 and —1 levels in the 
El and E2 symmetry species, correspond to K > and 



K < respectively. Because of nuclear-spin statistics, 
in the ground vibrational state the nondegenerate lev- 
els of J =cvcn, K — are only allowed to possess the 
overall symmetry Al, Bl, whereas levels with J =odd, 
K = are only allowed to possess the overall symmetry 
A2, B2. The K = doubly degenerate levels of El and 
E2 symmetry are denoted by +1 levels, i.e. by El + 1, 
E2 + 1 levels. The exact ordering of the different sym- 
metry levels within a certain J, K level is determined by 
the relative contributions of the h^ v and hi v parameters 
(see for example Fig. 3 of Ref. [7J. The internal motions 
are strongly coupled to the overall rotation resulting in 
a strong dependence of the torsional-wagging energies on 
the K quantum number. Thus the level ordering may 
differ from one if -ladder to another. This turns out to 
be important for obtaining large enhancement factors, as 
it may result in closely spaced energy levels with a differ- 
ent functional dependence on /1 which are connected by 
a symmetry allowed transition. 

The panel on the right hand side of Fig. [2] shows an 
enlarged view of the J — 2, K — and J = 1, K = 1 lev- 
els, with all symmetry allowed transitions assigned with 
roman numerals. Note that transitions with AJ — in 
the K = manifold are not allowed. The transitions 
labeled by iii,iv,vi,vn,vili and x are of particular inter- 
est as these connect the closely spaced levels of different 
if-manifolds and have an enhanced sensitivity to a vari- 
ation of fi. A similar enhancement occurs for transitions 
between the J = 5, K = 1, and J = 4, K = 2 levels as 
well as for transitions between the J = 6, K = 1, and 
J = 5, K = 2 levels. In what follows, we will outline the 
procedure to calculate the sensitivities of these transi- 
tions. The resulting sensitivity coefficients are presented 
in Table HJ and Table [En] and discussed in SecHVl 

In the present work, we use the group-theoretical high- 
barrier tunneling formalism developed for methylamine 
by Ohashi and Hougen [8], which is capable of repro- 
ducing observations of the rotational spectrum of the 
ground vibrational state of CH 3 NH 2 to within a few 
tens of kilohertz (9j [10] . The high-barrier formalism as- 
sumes that the molecule is confined to one of n equivalent 
equilibrium potential minima for many vibrations, but 
that it occasionally tunnels from one of these n minima 
to another. The formalism fits in between the infinite- 
barrier approximation, where no tunneling splittings are 
observed, and the low-barrier approximation, where the 
present formalism breaks down. A backward rotation of 
the whole molecule is introduced to cancel the angular 
momentum generated by one of the large amplitude mo- 
tions - the so-called internal axis method - requiring the 
usage of extended group ideas. The reader is referred to 
Refs. [7HTT] for a detailed description of the high-barrier 
tunneling formalism and the used Hamiltonian. 

Table Q] lists the molecular constants used in our cal- 
culations. It includes three types of parameters: 'non- 
tunneling' or pure rotational parameters; parameters as- 
sociated with pure methyl torsion motion (odd numerical 
subscripts n); and parameters associated with the NH 2 
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Figure 2. Energy of the lowest rotational levels in the ground state of methylamine ( 12 CHg 4 NH 2 ). The levels are denoted by 
J, K and the overall symmetry. The panel on the right hand side of the figure shows an enlarged view of the J = 2, K — and 
J — 1,K — 1 levels, with all symmetry allowed transitions assigned with roman numerals. The sensitivity of these transitions 
are listed in Table |ITT| The two transitions that are designated with bold arrows and are labeled by III and IV have sensitivity 
coefficients equal to — —19 and = 2, respectively. 



wagging motion (even numerical subscripts n). The ob- 
tained /i-scaling relations for the different parameters of 
the high-barrier tunneling formalism of methylamine are 
listed in the rightmost column of Table [TJ In the next 
sections, we will discuss the scaling relations for the low- 
est order parameters, the scaling relations for the higher 
order parameters, and the problems encountered in deter- 
mining these, are discussed in the supplementary online 
material |12) . 



III. SCALING RELATIONS OF THE 
MOLECULAR PARAMETERS 

We will use two different approaches for determining 
the //-dependence of the molecular constants that appear 
in the Hamiltonian: 

(i) The first approach is based on the fact that the 
tunneling model essentially assumes that for each large- 
amplitude tunneling motion the system-point travels 
along some path in coordinate space. In zeroth approxi- 
mation, we may represent each large amplitude motion as 



a one-dimensional mathematical problem after parame- 
terizing the potential along the path and the effective 
mass that moves along it. Thus, for each large ampli- 
tude motion, we will set up a Hamiltonian that contains 
one position coordinate and its momentum conjugate. 
The parameters of this one dimensional Hamiltonian may 
be connected with the observed splittings which are fit- 
ting parameters of the high-barrier tunneling formalism. 
The parameters of the one-dimensional Hamiltonians are 
functions of the moments of inertia and the potential bar- 
rier only, and their /i-dependence can be found in a sim- 
ilar fashion as was done for methanol and other internal 
rotors [31 0]. Application of this approach is straight- 
forward in the case of the leading tunneling parameters 
of methylamine but some ambiguities appear for the J 
and K dependences of the main terms, because there are 
several ways of representing these dependences in a one 
dimensional model. 

(ii) In the second approach, we use the spectroscopic 
data of different isotopologues of methylamine to esti- 
mate the dependence of the tunneling constants. In anal- 
ogy with methanol, we expect the tunneling splittings to 
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Table I. Molecular parameters, P s , of the ground torsional state of methylamine CH 3 NH 2 [9], and their sensitivity to a variation 
of the proton-to-electron mass ratio p, defined as K^ B = p-Qf 3 -- All molecular parameters are in MHz, except p and pK, which 
are dimensionless. 
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a These parameters do not involve tunneling motions. 

b These parameters arise from the NH 2 inversion tunneling motion. 

c These parameters arise from the CH 3 torsional tunneling motions. 



follow the formula [3]: 

WspUtting — , e ai ^ Ired . (2) 
V ^ red 

This formula originates from the semi-classical (Wentzel- 
Kramers-Brillouin (WKB)) approximation that assumes 
that the effective tunneling mass, represented by I re d, 
changes with isotopic substitution, but that the barrier 
between different wells remains unchanged. This ex- 
pression was successfully applied to the J = 0, K = 
A — E splittings and the J = 1, \K\ = 1 splittings in 
methanol [3]. In methylamine, the h nv parameters cor- 
respond to the splittings in the J = 0, K = due to 
tunneling between framework |1) and framework \n) (the 
set of frameworks represent the equivalent potential wells 
between which the system can tunnel), and application 
of the WKB approach to these parameters is straightfor- 
ward. Moreover, since in fact all tunneling parameters 
in methylamine may be related to the same type of over- 
lap integral as the h nv parameters, we may expect that 
the isotopologue dependence of all tunneling terms can 
be described by Eq. ([2]). Unfortunately, ambiguities ap- 
pear again when we apply this approach to higher order 
terms in the methylamine Hamiltonian. These ambigu- 
ities are connected to the fact that vibrational basis set 
functions \n) localized near various minima are not or- 
thogonal, but in fact have nonzero overlap integrals with 
each other. The correlation problems that arise in the 
high-barrier tunneling formalism due to nonorthogonal- 



ity of the basis functions are discussed in some detail in 
Ref. [IB]. The main consequence which affects the iso- 
topologue approach is that there may be 'leakage' from 
one parameter to another; each fitted parameter appears 
as a sum of the 'true' parameter value plus a small lin- 
ear combination of all other parameters with a coefficient 
that goes to zero when the overlap integral goes to zero. 
While this effect should be insignificant for the main tun- 
neling parameters of methylamine, it may be important 
for higher order terms because even a small 'leakage' of 
the low order parameters may be comparable in magni- 
tude with the 'true' values of the higher order parameter. 

In order to verify the mass dependence coefficients for 
the parameters of the methylamine Hamiltonian, we have 
refitted available data on the CH 3 ND 2 [H], CD 3 NH 2 [H] 
and CD 3 ND 2 |16j isotopologues of methylamine using the 
high-barrier tunneling formalism. Unfortunately, the 
amount of data available in the literature was rather lim- 
ited; 66 transitions for CH 3 ND 2 [2], 41 transition for 
CD 3 NH 2 [15] and 49 transitions for CD 3 ND 2 [16]. There- 
fore, many of the higher order terms were not determined 
in the fits, while some low order parameters were de- 
termined with a few significant digits only. As a re- 
sult, it was possible to obtain the /x-dependence of the 
main tunneling parameters h 2v and h 3v only. In order to 
obtain information on higher order terms, wc have un- 
dertaken a new investigation of the CH 3 ND 2 spectrum 
with the Kharkov millimeter wave spectrometer. The 
newly obtained dataset for CH 3 ND 2 contains 614 transi- 
tions, comparable to the number of microwave transitions 



5 



available for CH 3 NH 2 (696 transitions). The CH 3 NH 2 
and CH 3 ND 2 fits have an almost equal number of var- 
ied parameters and obtained similar weighted root-mean- 
square deviations. The results of the CH 3 ND 2 investiga- 
tion will be published elsewhere [17) , here we will use only 
those results necessary for obtaining the scaling relations. 



A. Pure Rotational Constants 

The pure rotational or 'nontunneling' parameters in 
the model are connected to the usual moments of inertia 
of the molecule and to the centrifugal distortion param- 
eters. Therefore, we will assume the same ^-dependence 
for these parameters as used for methanol [3]. 



B. 



CH 3 torsion and the hs v parameter 



The h 3v parameter in the high barrier-tunneling for- 
malism corresponds to a pure torsion motion. The quan- 
tity |3/i3„| may be related to the usual E-A internal ro- 
tation splitting in a molecule that contains a group of 
C 3v symmetry. Assuming that the potential barrier is 
described by a cosine function and taking the moment 
of inertia of the methyl top to represent the mass that 
tunnels, we may set up a one-dimension internal rotation 
Hamiltonian 



FjPl + -7T (1- cosn 7 ) 



(3) 



with n = 3 for a threefold barrier, p 7 = —id/d"/ is the 
angular momentum operator associated to the internal 
rotation coordinate, F 7 is the internal rotation parame- 
ter and V3 the barrier height. Using a value for F 1 de- 
rived from the molecular constants, we may fit the barrier 
height V3 to the observed value for |3/i3„| and estimate 
the /i-dependence of h 3v . 

In the used axis system, the off-diagonal contribution 
to the inertia tensor is represented by the si parameter. 
For methylamine, this parameter is set to zero as it is 
not required by the fit. Thus, we may assume that the 
methyl top axis coincides with the principal axis a, and 
p = I 7 /I a > and F 1 = C conv /((l - p)I 7 ), with C conv being 
a conversion factor {C conv — 16.8576291 amu A cm -1 ). 
Using values for p and F (recalculated from rotational 
parameters) from Table m we obtain J 7 = 3.18 amu A 2 
and F 7 = 15.12cm -1 (ab initio value 15.1684cm" 1 [Tg]V 
The value for J 7 is close to the expected one which sup- 
ports the validity of the present analysis. Now, using this 
value for F 7 and the value for h 3v from Table [ij a fit to 
Eq. |3| yields the effective barrier height V3 = 683.7 cm -1 
(ab initio value 708.64cm -1 [18 J. The one-dimensional 
model with this value for V3 predicts values for the first 
torsional band and the A — E splitting in the first ex- 
cited torsional state that are in a good agreement with 
the observed values (269 cm -1 versus 264 cm -1 [T5] for 
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Figure 3. hi v and /i3„ parameters as a function of the reduced 
moment of inertia for the torsional and inversion motion for 
four different isotopologues of methylamine. The solid lines 
are fits according to Eq. |2| through the values of CH 3 NH 2 
and CH 3 ND 2 . 



the band origin and 186 GHz versus 180 GHz [19] for the 
splitting in the v t =X). All this indicates that the one 
dimensional model is physically sound and sufficiently 
accurate for our purposes. 

Finally, we obtain the /^-dependence of h 3v via 



k: 



p d(h 3v ) 
h 3v dp 
F 1 d(h 3v ) 
h 3v dF 7 



(4) 



where we have used the fact that F 1 scales as /z _1 , i.e., we 
assume that the neutron mass has a similar variation as 
the proton mass. The numerical evaluation d(h 3v )/dF 7 
using Eq. Q yields K 1 ^- = -4.66. 

In the upper panel of Fig. |3j the value of the h 3v 
parameter is plotted as a function of the reduced mo- 
ment of inertia, Ij ed — C conv jF 1 for 4 different iso- 
topologues of methylamine. As mentioned, the quan- 
tity |3/i3„| corresponds to the usual A — E internal ro- 
tation splitting in a methyl top molecule, hence, we 
expect the tunneling splitting to follow Eq. ([2]). The 
solid line in the upper panel of Fig. [3| corresponds to 
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a = 10.3 THz (amu A 2 ) 1/2 and ai = 7.84 (amu A 2 ) 1/2 , 
obtained using the CH 3 NH 2 and CH 3 ND 2 data. The re- 
duced moment of inertia is directly proportional to p. 
Thus, the sensitivity coefficient is given by: 



K n 3v = 



a, d(h 3v ) 




(5) 



From the above expression, we find for the h 3v parameter 
of CH 3 NH 2 a sensitivity coefficient of K^ 3v = —4.64, in 
excellent agreement with the value found from the one- 
dimensional Hamiltonian model. 



C. Inversion and the h-2 V parameter 

The interpretation of the h 2v parameter in terms of 
an effective mass moving in a one-dimensional effec- 
tive potential is not straightforward. For instance, ab 
initio calculations of the kinetic parameter for the in- 
version motion in the equilibrium geometry range from 
9.6017cm" 1 [IH] to 26.7291cm" 1 [20], while the barrier 
height in different studies varies from 1686 cm" 1 [21] to 
2081 cm -1 [22 . Since the system needs to tunnel six times 
in order to return to its initial configuration, we will treat 
this large amplitude motion as a six-fold periodic well 
problem, following Ohashi et al. 23 . Furthermore, we 
assume that the potential along the path can be repre- 
sented by a rapidly converging Fourier series. Thus, we 
use Eq. ^ with 7 replaced by r and n = 6 as a zeroth 
order model. The effective inversion-torsion constant F T 
and barrier height V§ can be determined from the split- 
tings in the ground state and NH 2 wagging band origin 
(780cm" 1 [2U). From this, we obtain F^.I^ct^ 1 and 
Ve =2322 cm" 1 , close to the values obtained by Ohashi et 
al. [23]- Following the same procedure as for h 3v , we 
obtain the p dependence of h 2vi K^ 2v = —5.49. 

In the lower panel of Fig. |3j the value of the hiv pa- 
rameter is plotted as a function of the reduced moment 
of inertia, I^ ed = C conv /F T for 4 different isotopologues 
of methylaminc. The solid line in Fig. [3] corresponds to 

2 and ai = 7.35 (amu A 2 )~ 1/2 
data. From 

this fit, we find for the h 2v parameter of CH 3 NH 2 a sen- 
sitivity coefficient equal to K^ 2v — —5.48, again in ex- 
cellent agreement with the one-dimensional Hamiltonian 
model. 



a = 44.4 THz (amu A 2 ) 1 
obtained using the CH 3 NH 2 and CH 3 ND 2 



system by the two large amplitude motions. In methy- 
laminc, q 2 and r 2 represent the interaction of the angular 
momentum generated by the NH 2 inversion and the 'cor- 
rective' 7r/3 rotation of the CH 3 group with the J z and 
J y components of the total angular momentum, respec- 
tively. It can be shown in different ways that q 2 has the 
same dependence on p as h 2v . For instance, it follows 
from a study of the correlations between the q 2 , 93 and 
p parameters carried out by Ohashi and Hougen [5]. In 
methylamine, two possible choices exist for p. p can be 
chosen such that Coriolis coupling due to the inversion 
plus corrective rotation is eliminated (q 2 fixed to zero) , or 
such that Coriolis coupling due to the internal rotation 
of the CH 3 group is eliminated (93 fixed to zero). These 
two choices result in a difference Ap = (3/ir)q 2 /h 2v [8]. 
Since p is in both cases a (dimensionless) ratio between 
different moments of inertia and independent of p, the 
above equation implies that q 2 and h 2v should have the 
same /i-dependence. 

From the CH 3 ND 2 isotopologue data, a sensitivity co- 
efficient K q2 — —5.53 was found, in good agreement 
with the K q2 = —5.50 obtained from the one-dimensional 

model and close to the value for K^ 2 " . The r 2 term is 
expected to have the same /u-dependence as q 2 . We were 
not able to check the isotopologue dependence for this 
term, since it was not required by the CH 3 ND 2 fit. 



E. Higher order terms 

The /i-dependence of the higher order terms, including 
the J and K dependences of the h 2v and h 3v parameters, 
was determined in a similar fashion (see the online mate- 
rial to this paper [H]). Unfortunately, some ambiguities 
and discrepancies between the different approaches ap- 
peared in the determination of the scaling relations for 
some higher order terms, which is reflected by the rather 
large error for these parameters (see Sec. IV). This is 
not a serious concern as the higher order tunneling pa- 
rameters only marginally affect the coefficients of the 
considered transitions. 



IV. SENSITIVITY OF SELECTED 
TRANSITIONS 

Using the scaling relations for the high-barrier tunnel- 
ing Hamiltonian determined in the previous section, we 
are now able to calculate the sensitivity coefficient of any 
desired transition in the ground state of methylamine. In 
order to do numerical calculations, we rewrite Eq. (JlJ as 



D. qi and parameters 

The linear terms q and r correspond to the interaction 
of components of the total angular momentum with the 
angular momentum generated in the molecule-fixed axis 



TSVmr, van mn i a\ 

» = ~^2ev ' (6) 

with v mn the transition frequency between state m and 
n for the present value of p and u^ n the transition fre- 
quency when p is replaced by p(l ± e) with e a number 



7 



Table II. Transitions in methylamine (CH 3 NH 2 ) that are detected in astro-physical objects in our local galaxy as listed in Lovas 
et al. |25j . The fourth column lists the transition strength multiplied by the electric dipole moment, /i e , squared. The last 
column lists the sensitivity of the transitions to a possible variation of the proton-to-electron mass ratio. 
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much smaller than 1 (in our calculations, we typically 
use e = 0.0001). v mn is calculated using values for the 
molecular constants as listed in Table [ij v^ nn and 2/~„ are 
calculated using the molecular constants scaled accord- 
ing to the relations that were determined in the previous 
section. 

We have calculated the coefficients for all rotational 
transitions in the ground state of methylamine with J < 
30 and K a < 15 and v mn below 500 GHz. The two largest 
coefficients Ku w — 19 and Ku « 24 were found for the 
liA2 <- 2 Q Al and 13 3 £T + 1 +- 12 4 E1 + 1 transitions at 
2166 MHz and 1458 MHz, respectively. 

In Table [TTJ the transitions of methylamine that are 
detected in astrophysical objects in our local galaxy are 
listed together with their transition strengths and sen- 
sitivity coefficients. Table |III| lists transitions involving 
levels that have an excitation energy below 10 cm -1 , i.e., 
transitions involving levels that are expected to be pop- 
ulated in cold molecular clouds. The rotational transi- 
tions labeled with an asterix have recently been detected 
by Muller et al. [6] via absorption in a cold cloud at a 
redshift z = 0.89. Due to their rather large transition 
frequency their sensitivity coefficients are only slightly 
enhanced. The transitions in Table IHII that are labeled 
by the roman numerals i-x, correspond to transitions in 



the J = 1, K = 1 and J = 2, K — levels that are shown 
in the right hand side panel of Fig. [2] The transitions 
labeled by I and II, correspond to transitions between 
the levels of if -doublets, hence these have sensitivities of 
approximately —1. The transitions labeled by V and ix 
are transitions between levels which splittings are signifi- 
cantly affected by tunneling motions. The sensitivities of 
these transitions are on the order of —5, comparable to 
the sensitivity of the h,2v and hs v parameters. The tran- 
sitions labeled by in,iv,Vl,Vll,Vlll and X are of particular 
interest as these are transitions between levels that differ 
in overall rotational energy as well as torsional-wagging 
energy. Consequently, cancellation may take place that 
lead to an enhancement of the sensitivity coefficients. Of 
these, the transition labeled by in has the smallest tran- 
sition frequency (2166 MHz) and the highest sensitivity 
coefficient (if M = —19). The transition labeled by iv at 
4364 MHz has a sensitivity coefficient equal to = +2. 

The estimated uncertainties of the coefficients are 
quoted in brackets in units of the last digits. There are 
two sources of the uncertainty in the coefficients: (i) 
the uncertainty in the determination of the molecular 
constants and (ii) inexactness of the scaling relations of 
the Hamiltonian parameters including errors due to ne- 
glecting the fi dependence of the torsion-wagging poten- 



Table III. Transitions in methylamine (CH 3 NH 2 ) involving levels with an excitation energy lower than 10 cm" 1 . The fourth 
column lists the transition strength multiplied by the electric dipole moment, /i e , squared. The last column lists the sensitivity 
of the transitions to a possible variation of the proton-to-electron mass ratio. The transitions labeled by roman numerals 
correspond to the ones depicted in Fig. [2] The transitions labeled with an asterix have recently been detected by Muller et 
al. [H] in a cold cloud at z — 0.89. 
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tial of the molecule. We have assumed the error in the 
scaling coefficients to be ±0.02 for the rotational param- 
eters, ±0.1 for the tunneling parameters h 2v , h^ v , q 2l r 2 
and ±1 for higher order tunneling terms. Since the un- 
certainties for the measured transition frequencies in the 
ground torsional state of methylamine are less than 10 -4 
(and below 5x 10~ 6 for the low- J transitions of interest in 
the present study Q2]), we assume that the main errors 
in sensitivity coefficients are due to inexactness of the 
scaling relations of the Hamiltonian parameters. There- 
fore, similarly to the procedure adopted in Ref. [5], the 
coefficients were calculated taking either the upper or 
the lower bound for the scaling relations, corresponding 
to the upper and lower bounds of the assumed uncer- 
tainties. The difference was taken as an estimate of the 
uncertainty of the coefficients. In spite of the large 
uncertainties of the scaling relations for the higher order 
terms, the resulting errors in the coefficients of the 
different transitions are below 3%. To test the influence 
of the uncertainties in the scaling relations of the higher 
order terms, we have performed an additional calcula- 



tion where only the non-tunncling parameters and h 2v , 
h3 V , q 2 and r 2 were used to calculate the coefficients 
for different transitions. The difference between this cal- 
culation and the calculation with the full set of scaling 
relations was less than 1.7%, i.e., within the uncertainties 
presented in Table |IT| and |III| 

It is interesting to note that almost identical values for 
the sensitivity coefficients are obtained by using an equa- 
tion that directly connects the sensitivity coefficient of a 
transition with the sensitivity coefficients of the Hamil- 
tonian parameters: 



K Vn 



-UK 



'p. 



dE n 8E n 



dP„ dP s 



where 



dE m 
dP„ 



(m\d s \m) 



(7) 



(8) 



is the derivative of the energy level E m with respect to the 
Hamiltonian parameter P s used in the program to build 



9 



up the least-squares-fit matrix, and K^ s is the sensitivity 
coefficient with respect to the s-th Hamiltonian parame- 
ter. Eq. Q is based on the assumption that the energy of 
state \m) may be represented as E m — J2 S Ps( m \O s \m) . 
This assumption is valid when the Hamiltonian depends 
linearly on the parameters, i.e. that the Hamiltonian may 
be written as H = ^2 S P s O s . The high-barrier tunneling 
Hamiltonian used for methylamine depends nonlinearly 
on p, but as = 0, the transition sensitivity coeffi- 
cients calculated using Eq. agree well with the results 
obtained by using Eq. the w 0.4% difference is at- 
tributed to the pk term, which is also non-linear and 
which scaling coefficient is nonzero. 

From Eq. ([7| it is seen that contributions to K^™™ 
from different Terms in the Hamiltonian are proportional 
to the relative contributions of these terms to the tran- 
sition frequency. From this fact, it is obvious that the 
resulting sensitivity coefficients are mainly determined 
by the largest terms in the Hamiltonian and uncertain- 
ties in the scaling relations for the high order parameters 
do not significantly affect our results. 

Eq. |7| , illustrates that the largest enhancement is ob- 
tained for transitions that connect two near degenerate 
levels that have substantially different dependences on 
fi. The different dependence on /i is provided when the 
two levels contain non equal contributions from different 
types of motions in the molecule. In that case, a tran- 
sition 'converts' one superposition of rotation-torsion- 
wagging motion to another superposition of rotation- 
torsion- wagging motion. A significant enhancement is 
obtained when a 'cancellation' takes place, i.e. when two 
levels have nearly the same total energy due to quanti- 
tatively different contributions from various types of mo- 
tion in the molecule. 

From Eq. it is possible to obtain an upper limit for 
the sensitivity coefficient that we may hope to find in the 
ground vibrational state of methylamine. Considering 
the main, low order terms, the maximum splitting due 
to the tunneling motions, i.e., the maximum torsional- 
wagging energy difference between levels n and to, may 
be roughly taken to be 4(/i 2 „ + hs v ). Large enhancements 
of the sensitivity are expected for transitions that convert 
a considerable fraction of this energy into rotational en- 
ergy. Using Eq. (|7| and the values and sensitivities of the 
molecular parameters as listed in Table [TJ the maximum 
sensitivity that we may hope to find is 



K, = Kf ± — (ih 2v [K^ - Kf] 



4h 3v [K^-Kf]) 



-l±64800/z^ 



(9) 



with K™* = — 1 (i.e. that the of a rotational param- 
eters) and v nm the transition frequency in MHz. The 
sensitivities obtained from our numerical calculations are 
indeed found within these bounds. 



V. CONCLUSION 

Spectra of molecular hydrogen in highly redshifted ob- 
jects have been used to constrain a possible variation of 
the proton-electron mass ratio \x since the 1970s. How- 
ever, as the observed absorptions in H 2 correspond to 
transitions between different electronic states, these are 
rather insensitive to //; the sensitivity coefficients 
are in the range (—0.01,-1-0.05) [26-28 . For this rea- 
son even the highest quality H 2 absorption spectra in- 
volving over 90 lines, observed with the large dish Keck 
Telescope [29] and the Very Large Telescope [30] yield 
constraints \Afi/fi\ at the level of only 5 x 10~ 6 . 

The notion that specific molecules exhibit an enhanced 
sensitivity to /i variation is changing the paradigm for 
searching drifting constants on a cosmological time scales 
from the optical to the radio domain. The use of 
the NH 3 inversion transitions in the microwave range 
that have ^-coefficients of —4.2, has led to much 
tighter constraints on A/i//i 31-34J. It was recently 
pointed out that microwave transitions in the methanol 
molecule (CH3OH) have sensitivity coefficients in the 
range (—42, +53) [5]. In this paper, we showed that 
the sensitivity of microwave transitions in methylamine, 
CH 3 NH 2 , are in the range (-19, +24). 

Methylamine is particularly relevant as it was recently 
observed at z = 0.8859 in the intervening galaxy to- 
wards the quasar PKS 1830-211 [BJ. The sensitivity coef- 
ficients of the observed transitions at 78.135, 79.008, and 



89.956 GHz transitions were calculated to be K„ 



-0.87 



for the first two and K„ 



-1.4 for the third transition, 



respectively (see Table III). These three methylamine 
lines have a mean radial velocity of uch 3 nh 2 = —6.2+1.6 
km s^" 1 [BJ. With |AX M | = 0.563 and the uncertainty in- 
terval Av = 1.6 km s _1 , we obtain a preliminary estimate 
of Ap/p: 



Ap 




Aw 


M 




cAK^ 



< 9 x 10~ 6 



(10) 



where c is the speed of light. 

A tighter constraint on Ajit//z is obtained from the 
comparison of ^ch 3 nh 2 with the radial velocity of the 
methanol line at 60.531 GHz also detected at z = 0.8859; 
vch 3 oh = —5.3 ± 0.5 km s _1 [6j. According to Ref. [3], 



this transition has a sensitivity coefficient 



-7.4. 



In this case we have |AA^,| = 6.5 and Av = 0.9 ± 1.7 
km s _1 , which yields |A/i//t| < 10~ 6 . This estimate con- 
tains an unknown input due to possible non-co-spacial 
distribution of CH 3 OH and CH 3 NH 2 . More robust con- 
straints on A/i//i are derived from observations of lines 
of the same molecule. In this approach the low frequency 
transitions of CH 3 NH 2 at 2166 and 4364 MHz would be 
particularly attractive as the difference of their sensitiv- 
ity coefficients is AKn « 21. 
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